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ABSTRACT 

In this paper a set of analytic formulae are presented with which the partial derivatives 
of the flux obscuration function can be evaluated - for planetary transits and eclipsing 
binaries - under the assumption of quadratic limb darkening. The knowledge of these 
partial derivatives is crucial for many of the data modeling algorithms and estimates 
of the light curve variations directly from the changes in the orbital elements. These 
derivatives can also be utilized to speed up some of the fitting methods. A gain of ~ 8 
in computing time can be achieved in the implementation of the Levenberg-Marquardt 
algorithm, relative to using numerical derivatives. 
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1 INTRODUCTION 

In recent years, the discovery and further characterization of 
transiting extrasolar planets (TEPs) has provided unique in- 
formation about the nature of planetary systems. The anal- 
ysis of a planet which periodically eclipses its host star yields 
the physical radius, the inclination and the mass of the sys- 
tem in addition to the parameters which are gathered from 
the radial velocity m easurements. Since the discovery of the 
first s uch system (see lCharbonneau et a l. 2000; B rown et al.l 
2001), more than 40 other TEPs were discovered around 
other stars. Currently operating ground-based surveys are 
producing numerous new discoveries. The doubling period 
of the number of kno wn TEPs is below one year. Moreover , 
existing (CoRoT, see lBarge et"aT]|2008l : IXlonso et al.ll2008l ) 
an d planned space-born e instruments (e.g. Kepler Mission, 
sec iBorucki et aT]|2007h are expected to yield hundreds of 
new discoveries of such systems, even with planetary radii 
comparable to that of Earth. Also, subsequent observations 
of a given transiting system can provide some informa- 
tion on the variations in the timing of successive transits 
and the light curve shape. These detecti ons can be used 
to constrain other p l anetary companions (lAgol et alj 120051 : 
ISteffen fc Agoll2007l : iHolman fc Murravll2005l) o r co-orbital 
companions (Trojans. see lFord fc Holmanll2007l ). 

In order to optimize the precision and speed of TEP ob- 
servations, a careful analysis of the light curves is required. 
The basis of such light curve analysis is to find an adequate 
model of planetary obscuration, which causes a small de- 
crease in the stellar flux. Since both the stellar and the 
planetary body can be well modelled by a spheroidal shape, 
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the decrease in the stellar flux can be estimated from the 
full or partial overlap of two circles. The star itself has a 
non-negligible limb darkening which depends on both the 
stellar properties and the photometri c band and q uantified 
by a small set of coefficients (see e.g. IClareul2004l . for such 
tables of limb darkening constants). At present, the most 
widely used models f or t his problem have been given by 
iMandel fc Ago! (|2002| ) and lCimenej i|2006l ). lMandel fc Agoll 
( 20021 ) calculate closed form expressions for the flux decrease 
assuming non-linear or quadratic limb dark ening (qu a ntified 
by 4 or 2 coefficients, respectively) while ICimene 1 (|2006l ) 
gives an infinite series where the limb darkening can be 
taken into account up to arbitrary order. In most cases the 
quadratic case is adequate because the photometric preci- 
sion of typical data is not good enough for the higher order 
limb darkening models to make a difference. 

Most data modeling algorithms, including the well 
known nonlinear Levenberg-Marquardt fitting method (see 
iPress et al.lll992l ) utilize the partial derivatives of the model 
function with respect to the model parameters. The un- 
certainties in the model parameters can be well ch arac- 
terize d by the Fisher information matrix (see e.g. iFinnl 
1992), also requiring knowledge of the same partial deriva- 
tives. Therefore, in the case of planetary transits, the para- 
metric derivatives of the flux decrease function can be ex- 
tremely valuable. Moreover, partial derivatives can be used 
to construct a set of uncorrelated parameters of the light 
curve which is preferred by most of the parameter fitting 
algorithms (e.g. Levenberg-Marquardt, downhill simplex, 
Markov Chain Monte Carlo). Also, the analysis of the par- 
tial derivatives with respect to the limb darkening coeffi- 
cients themselves yields a combination of these with which 
consistent sanity checks can be done verifying the stellar at- 
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Table 1. Exclusion cases of different occultation geometries. To 
figure out the respective case for a certain value of p and 2, check 
the relations step by step: if the current relation is true, that case 
can be assigned to the given values of p and 2, otherwise go to 
the next step. The final case G is when there is no obscuration. 
Cases with subscripts can only occur if the radius of the planet 
is not smaller than 1/2. 
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mospheric properties in an independent way. Finally, these 
derivatives can be used to directly calculate how the varia- 
tions in the orbital parameters affect transit timing and the 
shape of the light curve. 

In this paper, we present the partial derivatives of the 
flux decrease function assuming a quadratic limb darkening 
law. In the next section, the formalism and the derivatives 
are presented. In Section 3, we apply these derivatives to 
construct a set of well behaved parameterizations of transit- 
ing light curves which yield an always finite and moderate 
correlation between the adjusted quantities in all important 
cases and limits. The correlations between the limb darken- 
ing coefficients are also discussed. The results are summa- 
rized in the last section. 



W 2 = 

w ' 



(4) 
(5) 



where W — 6 — 271 — 72. In equation ((2| the terms Fo, Fi, 
Fk, Fe, Fu and F2 are only functions of the occultation 
geometry, namely the relative planetary radius p = R p / 
and the normalized projected distance z between the center 
of the star and the center of the planet. In equation (J2J) the 
functions K(-), E(-) and II(-,-) denote the complete ellip- 
tic integrals of the first, second and third kind, respectively. 
The variation in the occultation geometry yields 12 distinct 
cases of obscuration, which are summarized in Table [1] This 
table is an exclusion table and should be interpreted as fol- 
lows. For a given value of (p, z), the first relation (in step 1) is 
checked. If it is true, the appropriate case can be assigned to 
the geometry, otherwise the next relation should be checked 
and so on. The different cases are denoted by bold capi- 
tals for planetary radii smaller than 1/2 and capitals with 
a subscript which can only occur if the radius of the planet 
is greater than or equal to 1/2 In practice it would barely 
occur for planetary companions for earlier types of stars but 
it might happen in the cases when a Jovian planet transits 
a later main sequence star (M dwarf). For the actually most 
common planet-like applications (p < 1/2), the 7 major 
cases (A,. . . ,G) are in the order of growing distance between 
the geometrical centers of the planet and the star. In equa- 
tion ([2]) the expressions for the terms Fi (i = 0, 1, K, E, II, 2), 
k and n can be found in tables lATI and lA"2l in Appendix lAl 
after the appropriate geometrical case has been obtained. 
We should note here that equation © is completely eq uiv- 
alent with the equation found in lMandel fe Agol (|2002T ). in 
the first line of the second paragraph in their Section 4 at 
page L173. However, this expansion of equation ([2} clearly 
separates the terms which depend only on the limb darken- 
ing coefficients (W m ) and the terms which depend only on 
the occultation geometry (Fi). 



2 PARAMETRIC DERIVATIVES OF THE 
FLUX DECREASE 

In this section the partial derivatives of the flux decrease are 
presented. The surface brightness of a star as a function of 
the normalized distance ^ r 1 assuming quadratic limb 
darkening is given by the equation 



/(r) = l- ]T 7™. (l-v/l 



(1) 



where th e constants 71 and 72 q uantify the limb darkening. 
Recalling iMandel fe Agol (|2002l ). the relative apparent flux 
of an eclipsed star with a quadratic limb darkening can be 
written as / = 1 — A/ (assuming a unity flux out of the 
transit) , where flux decrease A/ can be calculated using the 
equation 

A/ = W0F0 + W 2 F 2 + (2) 
Wi[Ft + F K K(k) + F E E(k) + F n n(n, fc)]. 

In this equation the quantities Wi (i = 0, 1, 2) are only func- 
tions of the limb darkening coefficients, namely 

6 — 671 — I272 



W = 



W 



(3) 



2.1 Partial derivatives with respect to the limb 
darkening coefficients 

Since in equation ((2)1 the only quantities which depend on 
the limb darkening coefficients are Wo, Wi and W2, the par- 
tial derivatives of A/ with respect to 7 m (m = 1,2) can 
easily be obtained, namely 



9A/ 



d A + ^F 2 + 
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where the appropriate derivatives dWi/d'y m are the follow- 
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2.2 Partial derivatives with respect to the 
geometric parameters 

In equation (J5|, the terms F; explicitly depend on the rela- 
tive planetary radius p and the normalized distance z. There 
is also an implicit dependence via the complete elliptic in- 
tegrals since their parameters k and n are also functions 
of p and z. The derivation of these partial derivatives are 
quite straightforward for the non-degenerate cases, i.e. for 
the cases B, Bg, D and F since all of the appearing func- 
tions in these domains are analytic. For the other cases, the 
partial derivatives can be calculated as the appropriate lim- 
its, namely 
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(13) 
(14) 
(15) 
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For Aq and G, all of the derivatives are obviously 0, except 
for the case (p = 1, z = 0) when the partial derivatives do 
not exist. 

Utilizing the parametric derivatives of the elliptic inte- 
grals (see Appendix [B}, the final form of the partial deriva- 
tives of A/ with respect to z and p is 



dAf 
dg 



W F , g + WiF lig + W 2 F 2 , g + 
(F A - + F E )k g 



WiK(fc) 
WiE(fc) 
+ 



■K,g 



Fe, 9 + 



k 

F K k g 

fc(l-fc 2 ) 



+ 



F n n„ 



2n(n - 1) 



(19) 



+ 



F n kk„ 
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, F E k g 
F n n s 



(fc 2 -n)(l-fc 2 ) 2(fc 2 -n)(n-l) 

where g denotes either p or z, the appropriate geometric 
parameter. The expressions for Fo, g , Fi, g , Fk, 3 , Fe, 3 , F 2 , g , 
kg and n g can be figured out for all cases using the tables 



in Appendix El namely Table [Cl] and Table IC2l 

We note here that the computation of these derivatives 
are even more simple and faster than the computation of 
equation (|2j since equation (I19|l lacks the complete elliptic 
integral of the third kind for which evaluation requires most 
of the computing time. 



All of these applications assume a transiting planetary sys- 
tem on a circular orbit with a given semimajor axis (rel- 
ative to the radius of the star) a/R*, an impact parame- 
ter, 6 = (a/R-t) cos i, the planetary companion has a fixed 
mean motion of n — 2n/P and the transit occurs at the 
instance E. The relative radius of the planet is denoted by 
p = Rp/R t . Therefore, the distance between the center of 
the stellar and planetary disk has a time dependence, 



'(*)= -£- 



^JL^ sin 2 [n{t- E)]+b 2 cos 2 [n(t- E)]. 



(20) 



From now on we assume that the semimajor axis is relatively 
large, i.e. the distance can be approximated by 



E) 2 + b 2 . 



(21) 



In the following parts of this section, we first calculate the 
correlations between the limb darkening coefficients, assum- 
ing the orbital parameters and the relative planetary radius 
to be known. In the second part of this section, we construct 
a set of adjusted parameters which always yields finite (i.e. 
definitely smaller than unity) correlations between them in 
the cases of non-grazing eclipses. This is relevant for studies 
of transiting planets since equation (|21[) yields a unity cor- 
relation between a/7?* and b in the limit of p — > with or 
without limb darkening for all impact parameters. In the last 
part of this section we present an analytical calculation how 
the uncertainties in the light curve parameters depend on 
the photometric passbands, assuming a Jupiter-sized planet 
orbiting a solar-type s tar. In all o f these cases we use the 
Fisher matrix method (|Finnlll992h to obtain the uncertain- 
ties and correlations of the fitted parameters. This method 
gives the covariance matrix as 



(5a m 5a n ) = (r ) f 
where 



= E 



9 m /(a,ti)5„/(a, U 



(22) 



(23) 



/(a, ti) is the observed flux at the instance U, a = (ai, a 2 ) = 
(71, 72) or a = (ai, a 2 , 03, 04) = (p, a/R*, b, E) is the vector 
of the adjusted parameters (depending on the actual appli- 
cation) and dm = d/da m . Since we are interested only in the 
correlations between the parameters and we can expect uni- 
form uncertainties in the measurements and uniform data 
sampling, therefore F can be multiplied by any arbitrary 
constant and the sum in equation ((23} can be replaced by 
the integral 



oc / d m f(&,t)d n f(&,t)&t. 



(24) 



Here ti < E - {na/R-,y x (l + p) and E + (no/ii*)" 1 (l + 
p) < t 2 , assuring that the ingress and egress are completely 
observed, independently from the impact parameter. 



3 APPLICATIONS 

In this section we present three simple applications which 
utilize the partial derivatives of the flux decrease function. 



3.1 Correlations between the limb darkening 
coefficients 

Now, we determine the correlations between the limb dark- 
ening coefficients when the adjusted parameters are a = 
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0.4 0.6 
Impact parameter 

Figure 1. The dependence of the correlation between the limb 
darkening coefficients, 0(71,72) as a function of the impact pa- 
rameter. Thin lines are for p = 0.01, moderately thick lines are 
for p = 0.1 and thick lines are for p = 0.2. The continuous, 
long dashed, short dashed, dotted and dotted-dashed lines repre- 
sents the cases when limb darkening coefficients are 71 = 72 = 0, 
71 = 72 = 0.1, 71 = 72 = 0.2, 71 = 72 = 0.3 and 71 = 72 = 0.4, 
respectively. 



(71,72). It is easy to show that this correlation would only 
depend on the impact parameter, the planetary radius and 
the two limb darkening parameters themselves while it does 
not depend on the mean motion, geometrical ratio of a/R+ 
and the transit center time E. We have obtained these cor- 
relations for very small (p = 0.01), average (p — 0.1) and 
large (p = 0.2) planetary radii assuming limb darkening co- 
efficients between 0.0 and 0.4, while the impact parameter 
was varied between and 1. We found that the correlation 
is always negative, relatively large, i.e. |C(7i,72)| > 0.93 
and it strongly depends on the impact parameter. For these 
certain values, the correlation is plotted as a function of b 
on Fig. [1] It can easily be seen that the smallest correla- 
tion is around b w 0.7 — 0.8, almost independent of the limb 
darkening and the radius of the planet. 

Let us now calculate the optimal linear combination 
of the limb darkening coefficients which can be adjusted to 
yield no correlation. Define the parameters ui and U2 as 



cos if sin ^1/71 
— sin ip cos ip J \ 72 



(25) 



For simplicity, let us denote the above orthogonal matrix by 
O = Omn- It can be shown that the covariance matrix of 
(uifUz) and that of (71,72) are related to each other as 



(SukSuf) = Ohm {SjmS'yn) O n t. 



(26) 



To make the matrix {SutSui) diagonal, the rotation param- 
eter ip should be in coincidence with the orientation of the 
eigenvectors of {6^ m 5^ n ), namely 

2 (071671) - (0721372) 

We have obtained the optimal values of the rotation pa- 
rameter p for very small (p = 0.01), average (p = 0.1) and 
large (p = 0.2) planetary radii assuming limb darkening co- 
efficients between 0.0 and 0.4, while the impact parameter 
was varied between and 1. We found that like above, this 




0.4 0.6 
Impact parameter 

Figure 2. The optimal rotation parameter (in degrees) to avoid 
correlations between the limb darkening parameters as a func- 
tion of the impact parameter (see text for definition and further 
details). Thin lines are for p = 0.01, moderately thick lines are 
for p = 0.1 and thick lines are for p = 0.2. The continuous, 
long dashed, short dashed, dotted and dotted-dashed lines repre- 
sents the cases when limb darkening coefficients are 71 = 72 = 0, 
7i = 72 = 0.1, 71 = 72 = 0.2, 71 = 72 = 0.3 and 71 = 72 = 0.4, 
respectively. 



parameter mostly depends on the impact parameter. The 
results are plotted on Fig. [5] The usefulness of such a plot is 
somewhat limited if we have no a priori knowledge from the 
limb darkening coefficients themselves since the correlation 
between them and therefore the optimal rotation parame- 
ter depends on the actual values of 71 and 72. However, in 
practice if we have a hint for the planetary radius and the 
impact parameter, this angle can be estimated within a few 
degrees since the correlation depends more strongly on b and 
p than 71 or 72. If we do not know any reasonable value for 
p or b before the fit, an angle of ip w 35° — 40° is a plausible 
selection in general. 



3.2 Correlations between the light curve 
parameters 

In this subsection we investigate the correlations between 
the light curve parameters utilizing the previously obtained 
partial derivatives and Fisher information matrix method. 
The classical formalism of adjusting parameters of a tran- 
siting system uses the same parameters as in equation (1201) . 
namely the ratio a/R+, the impact parameter b and the in- 
stance of the center of the transit E as well as the radius 
of the planet, p = R p /R*. Since the flux decrease depends 
directly on the radius p and indirectly on a/ 7?*, b and E 
(assuming a fixed limb darkening), the partial derivatives 
of the light curve fit) = 1 — A/(t) can be obtained by us- 
ing equation (|19[) and the chain rule. These derivatives can 
then be plugged into equation (|24[) while the adjusted set 
of parameters will be a = (p, b, a/R*, E). We have obtained 
the correlations between these variables and we have found 
that the correlation between b and a/7?* tends to unity as 
the radius of the planet is decreased. This correlation is plot- 
ted as a function of the impact parameter for planetary radii 
p = 0.01,p = 0.1 and p = 0.2 and for various limb darkening 
parameters on Fig. [3] It can clearly be seen that C(b, a/R*) 
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Figure 3. Correlation between the adjusted values of a/R* and 
the impact parameter b for various values of planetary radii and 
limb darkening coefficients: Thin lines are for p = 0.01, moder- 
ately thick lines are for p = 0.1 and thick lines are for p = 0.2. 
The cases for different limb darkening constants are almost indis- 
tinguishable. 

is almost independent from the actual limb darkening. From 
further analysis of the correlation, it turns out that for small 
impact parameters (b < 0.5), C(b, a/R*) can be well approx- 
imated by 1 — p 2 /2. Therefore, for very small planets, this 
correlation can be undesirable and most of the fitting meth- 
ods would distort the results. 

At this point we have checked the correlati ons between 
an al ternative parametrization proposed by iBakos et al.l 
In that work the light curve was parametrized by 
the equation 

* 2 W= fi-V(l-6 2 )(*--E) 2 + 6 2 > (28) 



where (,/R* = n(a / R+) / y/ 1 — b 2 . This parameter is re- 
lated to the duration of the transit, namely (£/_R*) -1 = 
Tkuration/2, where Td ura tion is the time between the instances 
when the center of the planet crosses the limb of the star. 
Since the above parametrization of z 2 is linear in b 2 , we have 
chosen b 2 instead of b as an independent parameter. We have 
found that utilizing the parameter set a' = (p,b 2 , £/ E) 
yields practically no correlation between b 2 and C,/R+ for 
non-grazing eclipses and increases only up to unity near 
6 > 1 — p. This correlation is plotted on Fig. [4] for various 
planetary radii and limb darkening coefficients. 

We should m ention here that the recent work of 
ICarter et all i|2008r i gives an exhaustive analysis of the un- 
certainties and correlations for various kind of transiting 
light curve characterizations. Their work focuses on the an- 
alytical calculations for light curves with no limb darkening 
and compares these results with numerical derivations for 
the limb darkened cases. 

3.3 Uncertainties of the light curve parameters 

In this subsection we calculate the dependence of uncer- 
tainties of the light curve parameters a/R±, E and b (see 
equation I20|l and the fractional planetary radius p for var- 
ious photometric passbands from near-ultraviolet to mid- 
infrared. It is known that the limb-darkening parameters 



decrease for longer wavelengths, therefore the transits them- 
selves become shallower and flattened. Using the Fisher 
information matrix method, as described earlier gives a 
simple and straightforward way to obtain these uncertain- 
ties. Assuming a solar-type star - i.e. with metallicity of 
[Fe/H] =0.00, surface gravity of logg* = 4.44 (CGS) and 
atmospheric temperature of T B a — 5780 K - we estimated 
these uncertainties for photometric passbands u , g' , r' , i', 
z' , J, H and K, when such a star is transited by a hy- 
potetical planet with the orbital parameters of a/R± = 10, 
P = 3.67 days and p = R p /Ri, = 0.1. The appropriate limb 
darkening coefficients for ea ch filter have been obtained us- 
ing the tables provided by IClaretJ (|2004h . The results are 
plotted on Fig. [5] for various impact parameters (b = 0.2, 
b = 0.5 and b = 0.8). During these estimations, the transits 
are assumed to be observed with a cadence of 10 seconds 
and with a photometric precision of Am = 1.4mmag. Note 
that this photometric precision is attainable by ~ 1— 1.2 m 
class telescopes for relativel y bright stars (z' ~ 9.5 mag) and 
using such cadence (see e.g. IWinn et al.ll2007h . 

It can clearly be seen from the plots of Fig. [S] that the 
uncertainties for all of the parameters decrease for longer 
wavelengths; moreover, this decrement can reach a factor of 
~ 10 (between the near-ultraviolet and mid-infrared bands) 
for the planetary radius. We note here that these analytical 
results agrees well with numerical estimations (Joshua N. 
Winn, personal communication). 



4 DISCUSSION AND SUMMARY 

In this paper the partial derivatives of the flux decrease func- 
tion of exoplanetary transits (or stellar binary eclipses) has 
been calculated assuming a quadratic limb darkening law. 
These derivatives can then be applied for various analyses 
from which we have demonstrated the correlation analysis of 
the limb darkening coefficients and two of the known tran- 
sit light curve parameterizations. The most time consuming 
part of the evaluation of the flux decrease (and its deriva- 
tives) is the computing of the complete elliptic integrals. 
Therefore, the calculation of the partial derivatives does not 
increase significantly the total computing time of both. The 
presented analytical analysis of light curve parametrization 
is extremely fast comparing to such an analysis based on 
Monte-Carlo methods: the integral in equation (|24[1 should 
be calculated only once instead of the evaluation of the x 2 
function O(10 4 ) timesfl The knowledge of these derivatives 
also allows straightforward calculations about how the vari- 
ations in the orbital elements affect the light curves and the 
timings of the successive transits. 

Here we note that such derivatives are also helpful in 
the implementation of the Levenberg-Marquardt algorithm. 
Since this method requires the partial derivatives of the func- 
tion to be adjusted, these must be evaluated either analyti- 
cally or numerically. The numerical evaluation of the partial 
derivatives requires the computation of the original function 
twice in all directions of the parameter space. Therefore for 
such a problem like transit light curve fitting, the numerical 

which is necessary to obtain a reliable a posteriori distribution 
of the parameters 
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Figure 4. Correlation between the adjusted values of Q/R* = 
n(a/R t )/\/l — b 2 and the square of the impact parameter, 6 2 for 
various values of planetary radii and limb darkening coefficients. 
Thin lines are for p = 0.01, moderately thick lines are for p = 0.1 
and thick lines are for p = 0.2. The continuous, long dashed, 
short dashed, dotted and dotted-dashed lines represents the cases 
when limb darkening coefficients are 71 = 72 = 0, 71 = 72 = 0.1, 
7i = 72 = 0.2, 71 = 72 = 0.3 and 71 = 72 = 0.4, respectively. 



approximation requires approximately 10 times more com- 
putation time (note that in practice the gain will be less, 
~ 8 due to other overheads resulted by the computation 
of increased number of coefficients and the gain will clearly 
depend on the used programming environment and its fea- 
tures). Moreover - because the derivatives of a transiting 
light curve function lack the Lipschitz property of continu- 
ity - the numerical approach can also be unstable at the 
points of contacts. 

The routines for calculating both the transit decrease 
function and its derivatives are availably in Fortran77 and 
C languages along with the codes used to calculate the cor- 
relations between the limb darkening parameters and the 
light curve parameters. 

The recent review of ISouthworthl l|2008h concludes that 
the determination of some of the light curve parameters, 
especially the radius of the planet can be sensitive to the 
applied limb darkening model and its parameters, yielding a 
possibility of systematic errors. Therefore, analytic descrip- 
tion of transit light curves for other different limb darkening 
laws would also be a point of interest. 
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Figure 5 . Uncertainties of the geometrical ratio a / i?* , the impact 
parameter b, the fractional planetary radius p = Rp/R* and the 
transit center time E for a hypothetical transiting planet with a 
radius of p = 0.1 orbiting its Sun-like host star on an orbit with 
a semimajor axis of a/R* = 10. The uncertainties are plotted 
as the function of the photometric passbands, for various impact 
parameters. See text for further details. 
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APPENDIX A: COEFFICIENTS FOR THE 
CALCULATION OF THE FLUX DECREASE 
FUNCTION 

In this section, the coefficients required to evaluate the flux 
decrease function - as it is given by equation Q - are sum- 
marized. The evaluation of the quantities Fq, Fi, Fk, Fe, 
Fn, F2, n and k is done in two steps. First, using Table |A"T1 
a set of auxiliary variables should be evaluated, all of them 
are a function only of p and z, i.e. do not depend on the 
limb darkening coefficients. Note that for a given geomet- 
rical case, not all of these quantities have to be calculated, 
only those that referred to in the appropriate row of Ta- 
ble IA2I (Moreover, it might happen that for some of these 
equations the value of p or z are out of the allowed domain 
if they are used in the cases when there is no need for them.) 
Second, using Table lA"2l the quantities Fi, k and n should 
by calculated by substituting the previously obtained values 
of the auxiliary variables. 



calculated, only those that are referred to in the appropriate 
row of Table IC2I Second, using Table [C2l the quantities Fi :S , 
kg and n g should by calculated (where g is either p or z), by 
substituting the previously obtained values of the auxiliary 
variables. 



APPENDIX B: PARAMETRIC DERIVATIVES 
OF THE COMPLETE ELLIPTIC INTEGRALS 

The derivatives of the complete elliptic integrals of the first 
and second kind are the following: 



dK(k) 

dk 
AE(k) 



E(k) K{k) 



k(l-k 2 ) k 
E(k)-K(k) 



(Bl) 



(B2) 



dk k 
If the complete elliptic integral of the third kind is defined 
with the sign convention as 



t/2 



U(n,k) 



dtp 



(1 — n sin 2 <p) yl — k 2 sin 2 tp 



(B3) 



then its partial derivatives are the following: 



dU(n,k) 
dn 



dH(n,k) 
dk 



2(fc 2 -n)(n- 1) 
(n 2 - fc 2 )n(n,fc) 



E(k) + 



{k 2 -n)K(k) 



+ 



+ 



k 



n — k 2 



n 

E{k) 



k 2 - 1 



+ U(n,k) 



(B4) 
(B5) 



Throughout this paper we are using the sign convenction for 
IT(n, k) as it is defined by equation (IB3|) . 



APPENDIX C: COEFFICIENTS FOR THE THE 
CALCULATION OF THE PARTIAL 
DERIVATIVES OF THE FLUX DECREASE 
FUNCTION 

In this section, the coefficients required to evaluate the par- 
tial derivatives of the flux decrease function are summarized, 
which are needed to evaluate equation (|19p . The evaluation 
of these coefficients are also done in two steps. First, using 
Table IC1I a set of auxiliary variables should be evaluated, 
all of them are a function only of p and z (i.e. do not de- 
pend on the limb darkening coefficients). Note that for a 
given geometrical case, not all of these quantities have to be 
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Table Al. Au xiliary quantities for th e calculation of the flux decrease. Note that the quantities a, b, ko and k\ are defined similarly as 
were defined by Mandcl & Agol (2002) and the former two should not be confused with the same notations for the semimajor axis and 
the impact parameter. 
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Table A2. Coefficients for the flux decrease function. 



Step 


Case 


^0 


Fi 


Fx 






F 2 




k 


n 


1 


A 


P 2 


|(1-(P') 3 ) 











i 








2 


Ac 


1 


3 

















3 


B 


P 2 


2 
3 


CiCik 


CiCie 




Ip 2 (p 2 - 


h2z 2 ) 


l 4=pz 
V o 


a 


1 


Bf 


P 2 


Ti 











|p 2 (p 2 - 


h2z 2 ) 






5 


B G 


Go 


2 
3 




CgCgb 


CgGgii 


G 2 




jl — a 


a-1 

a 


6 


C 


P 2 


1 

I 4 
3 


&(i -V) 











2p 




7 


C r jr 


l 

4 











32 








8 


Co 


Go 




?p 2 ) ^16p(2p 2 - 1) 





G 2 




i 

2p 




9 


D 


P 2 





CiCik 






Ip 2 (p 2 - 


h2z 2 ) 


/ 4pz 
V l"^ 


a 


10 


E 


P 2 


Ti 











|p 2 (p 2 - 


h2z 2 ) 






11 


F 


Go 





GgGqx 


CgCge 


GgGcti 


G 2 




ll — a 


a-1 
a 


12 


G 



























Table CI. Auxiliary quantities for the calculation of the partial derivatives of the flux decrease function. 

C IK , P = +2p(l + 2(p 2 - z 2 )) C IE ,p = 14p(l - a) - 2(p - z)(z 2 + 7p 2 - 4) 

C IK ,z = "2z(5 + 2(p 2 - z 2 )) C IE ,z = 2z(l - a) + 2(p - z)(z 2 + 7p 2 - 4) 

C G K, P = -2(p 2 (12p + 21z) + z(z 2 - 4) + 2p(5z 2 - 6)) C G E, V = 4z(-4 + 21p 2 + z 2 ) 

C G k,z = 2p(4 - 7p 2 - lOpz - 3z 2 ) Ggb.z = 4p(-4 + 7p 2 + 3z 2 ) 

Go, P = f (2fc ) G 2 ,p = £ (2t 2 fc - 4zpsinfc ) 

G ,z = I (-2sinfc ) G 2 ,z = £ (2zpfc - (p 2 + z 2 + 1) sinfc ) 



Table C2. Coefficients for partial derivatives of the flux decrease function. 
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